Effects of disorder on the vortex charge 
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We study the influence of disorder on the vortex charge, both due to random pinning of the 
vortices and due to scattering off non-magnetic impurities. In the case when there are no impurities 
present, but the vortices are randomly distributed, the effect is very small, except when two or 
more vortices are close by. When impurities are present, they have a noticeable effect on the vortex 
charge. This, together with the effect of temperature, changes appreciably the vortex charge. In 
the case of an attractive impurity potential the sign of the charge naturally changes. 

PACS numbers; 74.25.Qt, 74.72-h 



I. INTRODUCTION 

Some time agoi it was proposed that the vortex in- 
duced by an external magnetic field in a type-II super- 
conductor should be electrically charged. This effect was 
proposed to occur since the chemical potential is ex- 
pected to be larger in the vortex core than in the bulk 
of the superconductor. It is energetically favorable for 
the electrons to lower their energy through the conden- 
sation energy and, since the vortex core is interpreted 
as being a normal region, the electrons tend to move to 
the bulk leaving a charge deficiency close to the vortex 
line. It is the electrochemical potential (the sum of the 
chemical potential and electrostatic energy) that is con- 
stant in the superconductor. Due to the circulating cur- 
rents around the vortex line the electrostatic potential 
is needed to compensate the centrifugal force due to the 
circular motion^. 

We should note that the initial understanding that the 
vortex core region is populated by normal electrons has 
been questioned. In particular, in the case of clean super- 
conductors, where the mean free path is much larger than 
the coherence length, the localized states bound to the 
vortex core are the result of Andreev scattering^. The 
core states are coherent superpositions of particle and 
hole states and interpreted as being the result of con- 
structive interference of multiple Andreev scattering from 
the spatial variation of the order parameter. Also it was 
shown that the main contribution to the supercurrent is 
originated in these states. 

Soon after the proposal of Khomskii and Freimuth it 
was suggested that the effect could be tested experimen- 
tally due to the dipole field created at the surface of the 
superconductor''. It has been claimed that the charge of 
the vortex has been measured using NMR in high tem- 
perature superconductors^. 

Theoretical studies on the existence of the vortex 
charge were carried out subsequentlySiSiS. In partic- 
ular, a relation was established between the vortex core 



charge and the vortex bound states for an s-wave vorteji^. 
In particular, in the quantum limit the influence of the 
bound states is important. This regime is reached when 
T/Tc < l/(fci?Co)i where Tc is the critical temperature, 
kp is the Fermi momentum and Co = vf/^q is the coher- 
ence length, corresponding to a regime where the thermal 
width is smaller than the level spacing. In this regime 
the particle-hole asymmetry in the local density of states 
(LDOS) has been related to the vortex charge^. The 
asymmetry results from the different effect the supercur- 
rent around the vortex has on the particle and hole wave 
functionsiS. However, in d-wave superconductors the low 
lying states are not localizedy*. As shown recently, it is 
the winding of the phase around the vortex and not the 
detailed decreasing of the gap amplitude near the vortex 
core that is ultimately responsible for the supercurrent 
and the nature of the low-lying states (bound states in 
the s-wave case) and their effect on many propertiesi2*i^, 
such as the vortex charge. 

Taking into account the screening of the vortex charge, 
the Friedel oscillation in the charge profile was obtained 
for the case of a s-wave vortes^, showing that the charge 
is screened but prevails with a somewhat reduced value. 
Other pairing symmetries were also considered®' show- 
ing that the existence of the vortex charge is universal. 
In all these studies the vortex is charged positively (elec- 
tron defficiency). This positive charge has been argued to 
be the cause of the Hall anomaly where the Hall conduc- 
tance changes sign when entering the superconducting 
phasei^. 

In the high magnetic field regime, where a Landau level 
description is appropriate, the low lying states are coher- 
ent through the vortex lattice^^ both in the s-wave^S, and 
in the d-wave cases^^. The coherent nature of the states 
originates gapless superconductivity due to the center of 
mass motion of the Cooper pairs and the action of a high 
magnetic fieldi^*ii. In this regime the vortex charge has 
not yet been studied. 

In d-wave superconductors it is likely that other order- 
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ings compete with the superconductivity. In particular, 
antiferromagnetismi^, or d-density waveaiS;. The vortex 
structure and, in particular, the vortex charge, have also 
been studied when there is competition between the var- 
ious order parametersi2iS2i2i-22-23 general, a small 
region around the vortex will have a non-vanishing order 
parameter which affects the density of states and in par- 
ticular the vortex charge. It may change sign from an 
electron defficiency to an electron abundance at the vor- 
tex core. Also, increasing the temperature, in a regime 
of parameters where the competing order is absent, the 
positive vortex charge is recovered^i. 

Also, since most systems have impurities, these affect 
both the motion of the quasiparticles through scatter- 
ing and through pinning of the vortices. Therefore it is 
necessary to study the effect of the impurities. Such a 
study has been carried out in the case of the d-density 
wave state24. It is clear that, at the very least, the effect 
of the addition of the impurities is to locally change the 
chemical potential. 

In this work we consider several causes for disorder and 
their effect on the vortex charge. The effect of a charge 
inside a superconductor is in general screened. In usual 
superconductors of the s-wave type, the dominant con- 
tribution is due to Thomas-Fermi like screening, as in 
the normal state2^. Due to the presence of the gap the 
Fermi surface Friedel oscillations are suppressed. The 
exponential screening acts on the scale of the Thomas- 
Fermi length which in general is much smaller than the 
coherence length, and therefore any charge is very small. 
However, in d-wave superconductors or in s-wave super- 
conductors in the vicinity of a vortex, the Friedel oscilla- 
tions are important since in the first case there are gapless 
states and in the second case there are states of an essen- 
tially normal character in the vicinity of the vortex corei. 
These oscillations act on a scale which is comparable to 
the coherence length, which in type-II superconductors 
is small. Therefore in the vicinity of the vortex core the 
screening effect, even though noticeable, does not change 
qualitatively the effect of the charge depletion. This is 
shown in Fig. 1 of ref.^ where the charge oscillations near 
the core, even though depressed, are still visible and only 
a quantitative change is observed. Therefore, for simplic- 
ity, we will neglect in this work the effect of screening 
since the results will be qualitatively the same. Also, as 
we will show later, the effect of the impurities has a local 
nature and only affects significantly the physical quan- 
tities near the vortex core where screening has not fully 
acted. Only far from the vortex core the screening of 
the Thomas-Fermi type will strongly suppress the charge 
oscillations^. 



II. VORTICES AT LOW TO INTERMEDIATE 
FIELDS 



Gennes equations Titp = ei/j where ip^{r) = {u* (r) , v* (r)) 
and where the matrix Hamiltonian is given by 
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The sums are over nearest neighbors {S — ±x, ±y on the 
square lattice); A(r) is the vector potential associated 
with the uniform external magnetic field B = V x A, 
the operator sg is defined through its action on space 
dependent functions, ssu{r) = u{r + S), and the operator 
fis describes the symmetry of the order parameter. It is 
convenient to perform a singular gauge transformation 
to eliminate the phase of the off-diagonal term Q in the 
matrix Hamiltonian. We consider the unitary FT gauge 
transformation H — > U~^TiU , where2& 
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with (^^(r)-|-(/)B (r) = 0(r). The phase field 0(r) is decom- 
posed at each site of the two-dimensional lattice in two 
components 0a (r) and (/'_B(r) which are assigned respec- 
tively to a set of vortices A, positioned at {rf}i=ij 



,Na, 



and a set of vortices B, positioned at {rf },i=i_7Vg . The 
phase fields 4'^=a.b are defined through the equation 



V X V(/.^(r) = 2^z^J(r-rf) 



(5) 



where the sum runs only over the /i-type vortices. After 
carying out the gauge transformation Q the Hamilto- 
nian reads 
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The phase factors are given by?^ '^s(^) ^ Ir ' '^^ 

and S(p{r) — (/)A(r) — (j)B{^), where ?ik^ = mv^ = 
hV4>fi — -A is the superfluid momentum vector for the 
/i-supercurrent. Physically, the vortices A are only vis- 
ible to the particles and the vortices B are only visible 
to the holes. Each resulting /i-subsystem is then in an 
effective magnetic field 



Consider the lattice formulation of a superconductor 
in a magnetic field. Let us start from the Bogoliubov-de 



B^'ff = X = B - 0oz^ <5^(r - rf) (7) 
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where each vortex carries an effective quantum magnetic 
flux 00 • For the case of a regular vortex lattica^SiSi, these 
effective magnetic fields vanish simultaneously on average 
if the magnetic unit cell contains two vortices, one of each 
type. More generally, in the absence of spatial symme- 
tries, as it is the case for disordered systems, these effec- 
tive magnetic fields B'^^^'^ vanish if the numbers of vor- 
tices of the two types A and B are equal, i.e. Na — Nb, 
and their sum equals the number of elementary quan- 
tum fluxes of the external magnetic field penetrating the 
system. 

The /i-superfluid wave vector k^(r) characterizes the 
supercurrents induced by the /.t-vortices. This vector can 
be calculated for an arbitrary configuration of vortices^i 
like 

i—l 

As we take the London limit, which is valid for low mag- 
netic field and over most of the H ~ T phase diagram 
in extreme type-II superconductors such as cuprates, we 
assume that the size of the vortex core is negligible and 
place each vortex core at the center of a plaquette. 

For the conventional s-wave case the operator charac- 
terizing the symmetry of the order parameter is constant 
fjs ~ J and the off-diagonal terms of the Hamiltonian © 
are then considerably simplified 

V s J 

(9) 

Note that in this case the phase of the off-diagonal term 
is eliminated. 

For the unconventional d-wave case the operator fis 
takes the form fjs — (— 1)'^"S5. With these definitions the 
d-wave Hamiltonian can be derived from the Hamiltonian 
© and reads 

V ^ s J 

(10) 

where the phase factor .45 (r) has the form 

As{r) ^ \ !^ {V4>A - V0b) • d\ 




In the Hamiltonian (|10|l and in Eq. ^Jthe vector 

a, = i(k^-kf) (12) 

acts as an internal gauge field independent of the external 
magnetic field'^^. The associated internal magnetic field 



b = Vxa^ consists of opposite A—B spikes fluxes carying 
each one half of the magnetic quantum flux (pQ, centered 
in the vortex cores and vanishing on average since the 
numbers of A- and _B-type vortices are the same. 

The solution of the BdG equations gives the spectrum 
and the wave functions of the quasiparticles. As the ef- 
fective magnetic fields experienced by the particles and 
the holes vanish on average, within the gauge transfor- 
mation we are allowed to use periodic boundary condi- 
tions on the square lattice (\E'(x -I- nL, y + mL) = '^{x, y) 
with n,m € Z) . The L x L original lattice becomes then 
a magnetic supercell where the impurities are placed at 
random and where the vortices are placed in such a way 
as to minimize their total energy. The disorder induced 
by the impurities in the system is then established over 
a length L. Thus in order to compute the eigenvalues 
and eigenvectors of the Hamiltonian we seek for eigen- 
solutions in the Bloch form ^ni^{r) = e~''^"{U*y^,V*^J 
where k is a point of the Brillouin zone. We diagonalize 
then the Hamiltonian e~^^'^TLe'^^'^ for a large number of 
points k in the Brillouin zone and for many different re- 
alizations (around 100) of the random impurity positions 
and of the correlated vortex positions. 

The density of states (DOS) and the local density of 
states (LDOS) for the cases of vortex disorder— and 
the combined effects of vortex disorder and impurity 
scattering2S: were studied recently. The disorder in gen- 
eral increases the DOS at low energies, by filling the gap 
in the s-wave case and originating a finite density of states 
at zero energy in the d-wave case. An approximate scal- 
ing regime was obtained in the last case. Also, the results 
for the LDOS are in qualitative agreement with STM ex- 
periments if most vortices are pinned at the impurity 
locations. 

The electron density is calculated in the usual way 
n{r) ^2Y,{\u^m^I{E.) + W^^l I{E.)) (13) 

i 

where Ei are the energy eigenvalues and f{Ei) is the 
Fermi function. 



A. Vortex lattice 

We consider first, for completeness, a regular distri- 
bution of vortices both for s-wave and d-wave pairings. 
At not very high fields, the vortices are sufficiently apart 
and the vortex charge is depleted from the vortex cores, 
as mentioned above. Therefore the vortices are positively 
charged. In Fig. ^ we consider the s-wave pairing case 
and in Fig. |2] we consider the d-wave pairing case. In 
this last case a checkerboard modulation of the charge 
density is seen. Both sets of results are presented for 
the case A ^ t. Linearizing the BdG equations defined 
in a continuum close to the nodes leads naturally to the 
definition of two velocities the Fermi velocity, vp, and 
a velocity va = Ag/pp, where pp is the Fermi momen- 
tum. This velocity denotes the slope of the gap at the 
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FIG. 1: (color online) Electron density n(r) — ribuik for s-wave 
symmetry in a regular vortex lattice. Here A — t, fi — —2.2t 
and ribuik = 0.377. 




FIG. 2: (color online) Electron density n(r)— ribuik for d-wave 
symmetry in a regular vortex lattice. Here A — t, fi — —2.2t 
and ribuik = 0.509. 



nod o^^i?? . Defining the anisotropy of the Dirac cone as 
the ratio ao = vf/va which in the lattice case translates 
to ajj — t/ Aq, the case considered above A — t is called 
the isotropic case. In many d-wave superconductors the 
anisotropy is actually stronger. We present for com- 
parison in Fig. 13 results for d-wave pairing considering 
A — 0.25i (note that in high- Tc materials the anisotropy 
is actually stronger of the order of t = 15 A). As the 
anisotropy increases the charge depletion decreases in 
depth but extends in area. 
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FIG. 3: (color online) Electron density n(r) — ribuik for d- 
wave symmetry in a regular vortex lattice. Here A — 0.25t, 
H — —2.2t and Wbuik = 0.36. 




FIG. 4: (color online) Electron density n(r) — ribuik for s-wave 
symmetry in a regular vortex lattice. Here A = /i = ~Q.lt 
and Wbuik = 0.965. 



We have considered a low density regime. We can vary 
the chemical potential (and the band-filling). In Fig. 0| 
we present the electron density for a situation close to 
half-filling for the s-wave case. In Fig. El we consider a 
case larger than half-filling. In this case there is a charge 
accumulation (or hole depletion). Indeed as explained 
before^ the density of the dominant carriers is depleted 
near the vortex core. Also, we have checked explicitly 
that at half-filling the vortex charge vanishes^ and the 
electron density is uniform. The results for the d-wave 




FIG. 5: (color online) Electron density n(r) — ribuik for s-wave 
symmetry in a regular vortex lattice. Here A = /i = 0.5t 
and Tibuik = 1.17. 



case show the same trend. 



B. Vortex disorder 

We consider now the case when the vortices are dis- 
tributed randomly due to some strong pinning effects, 
but neglect the effect of impurities on the motion of the 
quasiparticles. Therefore, the superfluid velocities are 
determined by a random distribution of the locations of 
the vortices, assumed static. 

In Fig. 1^ we consider the case of s-wave pairing and 
in Fig. d we consider d-wave pairing for a particular 
random distribution of the vortices (note that there is no 
sum over random configurations, like in the calculation 
of the LDOS). When a vortex is isolated the distortion 
of the charge profile is very similar to the lattice case. 
However, if two vortices are pinned nearby, the charge 
profile is significantly changed. In conjunction with the 
charge depletion, there is a sign reversal of the electron 
density. This is particularly visible in the s-wave case 
but also occurs for the d-wave symmetry. In the region 
where the two vortex cores are located there is a charge 
accumulation which leads to a local negative charge with 
respect to the bulk value. These fluctuations are of a 
similar order of magnitude as the charge depletion at the 
single vortices. However, the integration of the charge 
density in the neighborhood of the two vortices is still 
positive, despite the strong negative oscillations. The 
total charge in the case when two vortices are close by, 
even though still positive, is considerably smaller. 
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FIG. 6; (color online) Electron density n(r) — ribuik for s-wave 
symmetry in a disordered vortex lattice. Here the magnetic 
field is _B = 1/160 which corresponds to 10 vortices in a 40x40 
unit cell. Here A = t, = —2.2t and ribuik ~ 0.375. 




FIG. 7: (color online) Electron density n(r) —ribuik for d-wave 
symmetry in a disordered vortex lattice. Here the magnetic 
field is _B = 1/160 which corresponds to 10 vortices in a 40x40 
unit cell. Here A = t, /j. — —2.2t and ribuik ~ 0.507. 

C. Effect of impurities 

We consider now the effect of impurities. These are in- 
troduced at the Hamiltonian level, with the substitution 

Here U{r) is the potential due to the impurities placed 
randomly in the system. We model the disorder using 
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FIG. 8: (color online) Electron density n(r) — nbuik for d- 
wave symmetry with impurities and no vortices. Here A — t, 
H = -2.2t, U ^5t and ntuik = 0.509. 




FIG. 9: (color online) Electron density n(r) — ribuik for d-wave 
symmetry in a disordered vortex lattice with impurities. Here 
A = i, At = -2.2t, U = 5t, B ^ 1/160 and Wbuik = 0.505. 



the binary alloy modelS^. At each impurity site it costs 
an energy U to place an electron (it acts as a local shift 
on the chemical potential). The impurities are randomly 
distributed over a L x L periodic two-dimensional lattice 
and play the role of pinning centers for the vortices. In 
general, it is favorable that a vortex is located in the 
vicinity of an impurity^^. However, in this work we will 
be considering random distributions of the vortices and 



FIG. 10: (color online) Difference in the electron density 
between the cases with and without vortices shown in Figs. |H1 
andU 




FIG. 11: (color online) Electron density n(r) — ribuik for d- 
wave symmetry with impurities and no vortices. Here A = t, 
fi = —2.2t, U — —5t and ribuik ~ 0.513. 



impurities and study the electron density throughout the 
system with an arbitrary distribution of the vortices and 
impurities. Due to the non-homogeneous nature of the 
order parameter the BdG equations have to be solved 
self-consistently. 

The effect of the impurities is the expected one. If 
the potential is repulsive the electron density is lowered 
and if the potential is attractive the electron density is 
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FIG. 12: (color online) Electron density n(r) — ribuik for d- 
wave symmetry in a disordered vortex lattice with impurities. 
Here A = t, fj, ^ -2.2t, U = -5t, B = 1/160 and ribuik = 
0.509. 



To compare the effects of the impurities and the vor- 
tices, we show in Figs. |Hland|n|the electron density for 
U ^ 5t without and with vortices for a specific distribu- 
tion of the impurities. In Fig. ^|we plot the difference 
between the two cases. We see that even though the 
impurities have a strong (local) influence on the charge 
distribution, the contributions from the impurities and 
the vortices add up and the difference is not negligible. 
The same is observed for the attractive case U = —5t 
shown in Figs. 1111 1121 and 1131 Also, note that increas- 
ing the impurity potential for values \U\ > 5t does not 
change qualitatively the charge since |J7| = 5t is already 
a large value. 

Clearly, when the vortices are diluted and neglecting 
the effect of the impurities, the interaction between the 
vortices is not so important and the behavior of the sys- 
tem is not very different from an isolated vortex. A pos- 
sible exception is the regime of very high magnetic fields. 
But for low to intermediate fields, as seen from Figs. 
and 12 the distortion of the electron density occurs close 
to the vortex locations and the behavior is characteristic 
of a single vortex. The vortex charge has been studied for 
a single vortex, as mentioned above, considering a vortex 
with a quantum of flux. It is therefore also interesting to 
consider the influence of an impurity in the single vortex 
case. 




FIG. 13: (color online) Difference in the electron density 
between the cases with and without vortices shown in Figs, 
inland Hal 



increased leading to an electron accumulation and a sign 
reversal. There are sharp peaks at the impurity locations 
that mask the effect of the vortex charge in the material. 
If a vortex is pinned at an impurity site (as in most sys- 
tems they are) the effect of the impurity potential is quite 
strong. For moderate values of the potential the effect is 
smaller but still noticeable. 



III. SINGLE VORTEX 

We consider then the effect of an impurity in a s-wave 
vortex. We solve the Bogoliubov-de Gennes equations on 
a continuum-"^*^ introducing an impurity as a disc of small 
radius d centered at the vortex location. 

The general solution is obtained solving numerically 
the BdG equations. In the s-wave case these can be writ- 
ten as 

{He + U{r))u{r) + A{r)v{r) = eu(r) 
-(H: + U{r))v{r) + A*{r)u{r) = ev{r) (14) 

Here 

\ f H p \ ^ 
H,{r) = _ -V - -A + U{v) - Ef 
2m \i c ) 

where m is the electron mass, A is the vector potential, 
Ep the Fermi energy and V(y) is the potential originated 
from the impurity. 

The BdG equations are solved using the 
decompositions^^ : 

v\pM = Y.d\,,^'''^'i>,.Ap) (15) 

and choosing a gauge such that A.{p,Lp) = A{p)e~^"^'^ . 
appropriate for a vortex containing m flux quanta. Here 
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FIG. 14: (color online) Electron density of a single vortex as a 
function of distance for different values of T. The parameters 
are chosen as Ef = 1,V = 1.6,7? = 80. The energies are 
in units of the Fermi energy and the lengths in units of the 
inverse of the Fermi momentum. 

PL is an integer, j is the number of the zero of the Bessel 
functions in a disc of radius R and the normalized func- 
tions constitute a complete set over the zeros and 
are defined by 




Here is the Bessel function of order ^ and a^^^ is the 
j*'' zero of the Bessel function J^. The functions by 
construction, are zero at the border of the disc p = R. 
In the basis of the functions (jjj^^ we have to solve a set 
of equations for the c and d coefficients of the form 

Here the various components are diagonal in the angu- 
lar momentum T.f ■, = 5 a u'T^.^ and the same for 
the off-diagonal terms. In the case of strongly type-II su- 
perconductors and if we are interested in the low energy 
states, which are particularly relevant at small distances 
from the vortex core, the vector potential may be ne- 
glected. The various terms are then given by 

Jo 

and 

= / dpp(j)j^^{p)A{p)(j)j'^f,+^{p) (17) 
Jo 



FIG. 15: Vortex with a double flux 2$o. Left panel: A as 
a function of distance. Right panel: electron density as a 
function of distance 

where we have used that 

u{p)^ue{d-p) 

The gap function is obtained in the usual way 
Mp)^V u^Ap)v^Ap){l-'^f{E^,,)) (18) 

lJ.,i{\Ei\<uJo) 

with Ufj,Ap,(p) = e'■'^^'Uf_,,i{p) and w^,j(p, <p) = 
e^'^^^^^v^^p) where V is the attractive effective inter- 
action between the electrons. The BdG equations are 
solved self-consistently, as mentioned above. 

Considering first the case of a vortex enclosing a quan- 
tum of flux (to = 1) we show, for completeness, in Fig. 
Elthe charge profile close to the vortex core, considering 
first U — Q (no impurity Close to the vortex core the 
electron density decreases with respect to the bulk value. 
The effect is specially evident at low T. As the tem- 
perature increases the charge depletion decreases. This 
result is due to the Kramer and Pesch effect: the vor- 
tex core size increases as the temperature increases and, 
since the charge depletion may be related to the varia- 
tion of the gap function'^"', this derivative decreases as T 
increases and therefore the charge depletion decreases as 
T increases. As proposed in refi^ the electrostatic po- 
tential has contributions that are due to the difference of 
the gap function at the vortex with respect to the bulk 
value, due to the derivatives of the amplitude and of the 
phase of the gap function. Through Poisson's equation 
these dependences carry to the electron density. Note 
that at low T there are oscillations in the electron charge 
density. Also note that we are considering here the full 
quantum limit where we have access to the vortex struc- 
ture inside the vortex core. In the previous sections the 
method neglected the vortex core and we had no access 
to the true vortex core. 

As we saw in section II. B, when two vortices are close 
by there are also oscillations in the charge density that 
result, in this case, from the vicinity of two vortices. On 
a large scale (where the vortex core is averaged out) two 
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FIG. 16: (color online) Electron density and A as a function 
of distance for different values of U. The disc size is d = 0.96. 

vortices nearby may appear similar to a vortex contain- 
ing two flux quanta (note however that the energy of 
a double flux vortex is higher than two single quantum 
flux vortices). In Fig. E|we consider the case of a vortex 
containing two flux quanta (m = 2). The inner structure 
of the vortex is somewhat different^''^^ . The gap func- 
tion has a node near the origin and there are opposing 
currents in the same regime, as also obtained using the 
Andreev Hamiltonian. In the right panel of Fig. E| we 
show the electron density profile. As one approaches the 
core the charge is depleted, but very close to the location 
of the vortex the charge approaches the bulk value and 
changes sign. The trend is similar to the case studied in 
section II. B even though the effect is more pronounced 
in this last case. 

In Fig. El we consider the effect of the impurity po- 
tential at a low value of T. For positive values of C/ as ?7 
increases the charge is further depleted decreasing both 
the value of the electron density (it vanishes at the origin 
for a sufficiently large value U — 5) and extending the 
regime where the charge is depleted. A value of J7 = 5 is 
again similar to any larger value. If the impurity poten- 
tial is negative {U < 0) even small values of U like, for 
instance, U = —0.5 have a considerable effect on the vor- 



tex charge. The attractive potential accumulates charge 
at the vortex core. The effect of the impurity potential 
on the gap function if [/ > is also clear. Increasing 
U is similar to the effect of increasing the temperature 
(Kramer and Pesch effect). However, an attractive po- 
tential has a more profound effect. Even for a small value 
U — —0.5 a node appears in a way similar to the node of 
the vortex with two flux quanta (m = 2) previously con- 
sidered. As U decreases further for instance for U — —5, 
two nodes appear in the gap function. Therefore, there 
is a similarity between the attractive impurity case, the 
vicinity of two vortices and the multiple-flux vortex. The 
quantitative effect on the vortex charge is however dif- 
ferent, since the impurity potential is more effective in 
changing the signal of the vortex charge. 

IV. CONCLUSIONS 

Earlier treatments predicted an universal charge deple- 
tion at the vortex cores. Taking into account the com- 
petition in the d-wave case with other order parameters 
it has been determined that in some circumstances the 
vortices may be negatively charged (charge accumulation 
with respect to the bulk value). 

In this work we did not consider the effect ot other 
orderings but considered the inffuence of disorder. We 
focused on the effects of positional disorder of the vor- 
tices and on the effect of impurities. When two vortices 
are close by we found that strong fluctuations appear in 
the shared region of the vortices, that induce a smaller 
charge accumulation. Also, the addition of impurities 
changes the charge proflles. A small to moderate attrac- 
tive potential also changes the signal of the vortex charge, 
since it renormalizes locally the chemical potential in a 
straightforward way. 

The case of a vortex lattice in a very high mag- 
netic field, where the quasiparticles propagate coherently 
throughout the system and gapless superconductivity oc- 
curs, is a qualitative different state. In this regime a Lan- 
dau level description is adequate and a different behavior 
is found for many physical properties'^. Preliminary re- 
sults seem to indicate that the effect of the coherence on 
the vortex charge is to change the signal of the vortex 
charge: close to the vortices there is a charge accumula- 
tion instead of a charge depletion^. These results would 
then be in disagreement with possible explanations of the 
Hall anomaly as due to the positively charged vortices. 
Indeed the Hall anomaly is detected close to the normal 
phase, where in strongly type-II superconductors, like 
high-Tc materials, it is predicted that a Landau descrip- 
tion should be appropriateii. 
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